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We investigate the lowest state of a rotating, weakly-interacting Bose-Einstein condensate trapped 
in a harmonic confining potential that is driven by an infinitesimally asymmetric perturbation. 
Although in an ax;ially-symmetric confining potential the gas has an axially-symmetric single-particle 
density distribution, we show that in the presence of the small asymmetric perturbation its lowest 
state is the one given by the mean-field approximation, which is a broken-symmetric state. We 
also estimate the rate of relaxation of angular momentum when the gas is no longer driven by the 
asymmetric perturbation and identify two regimes of "slow" and "fast" relaxation. States of certain 
symmetry are found to be more robust. 

PACS numbers: 03.75.Fi, 05.30.Jp, 67.40.Db, 67.40.Vs 



I. INTRODUCTION 

A Bose-Einstein condensate of trapped alkali-metal 
atoms 1^] is an interesting system for studying phenom- 
ena connected with superfluidity. The behavior of this 
system under rotation has been studied extensively. Ex- 
perimentally vortex states in a two-component Bose- 
Einstein condensate were observed by Matthews et al. [|| . 
Also, Madison et al. |^ and more recently Abo-Shaeer et 
al. Q and Haljan et al. ^ have created vortices in a 
single-component condensate. 

Theoretical studies of this problem have been per- 
formed in the limit of weak interactions between the 
atoms |6 -16 , as well as in the Thomas- Fermi limit of 
strong interactions [T7|-|2^ . Up to now most experiments 
are in the Thomas-Fermi regime; however in the recent 
experiment of Gorlitz et al. in a cigar-shaped trap, 
the limit of weak interactions was reached transversely to 
the long axis of the trap. In addition, in the experiment 
of Ref. [Q, more than 100 vortices were created and as 
argued by Ho (l^, in such a system one can easily get 
to the regime of weak interactions. In the present article 
we restrict ourselves to the limit of weak interactions. 

As discussed by many authors [^-^ in this limit there 
is a degeneracy due to the harmonic confining poten- 
tial, which corresponds to the many different ways of 
distributing L units of angular momentum among N 
atoms. The goal, therefore, is to incorporate the in- 
teractions which lift the degeneracy and determine the 
lowest-energy state of the system. 

A basic question related to these systems concerns the 
lowest state for a given L and TV in a symmetric confining 
potential. References P-]!^] have dealt with this issue. 
Having understood the question of a symmetric confin- 
ing potential, one can then examine the ground state in 
the presence of a realistic trap, which has a small but 
finite asymmetry, or even in a more deformed trap. Our 
main goal is to investigate (i) the single-particle density 
distribution of the condensate when it equilibrates in the 
rotating frame, and (ii) the characteristic timescale for 



relaxation of angular momentum of the condensate when 
the trap stops to rotate and the gas is no longer driven 
externally. 

Our article is organized in the following way: 
We describe our model in Sec. II, examining how a 
small symmetry-breaking harmonic potential affects the 
lowest-energy state of a repulsive weakly-interacting 
Bose-Einstein condensate which is confined in a symmet- 
ric harmonic potential. In Sec. Ill we demonstrate that 
the state that minimizes the energy in the presence of a 
small asymmetric perturbation is the one given by the 
mean-field approximation, which is a broken-symmetric 
state. Finally, in Sec. IV we estimate the rate of relax- 
ation of angular momentum when the asymmetric drive 
stops to rotate, and in Sec. V we summarize our results. 



II. MODEL 



A. Symmetric confining potential 



Let us start with the Hamiltonian H, given by 



H 



E 



2M 



(1) 



with M being the atomic mass. In what follows we as- 
sume that the cloud rotates around the z axis, along 
which it remains in its ground state, so effectively we ex- 
amine a two-dimensional problem. Let us assume for the 
moment that the trapping potential is that of an isotropic 
harmonic oscillator of frequency uj±_ in the x-y plane. For 
this potential the single-particle energies e„ m are given 
by 



= (2rtr + |?7l| + l)TlLJ^^ 



(2) 



where is the radial quantum number, and m is the 
quantum number corresponding to the angular momen- 
tum. Let us also denote the corresponding single-particle 
states as ^„^,„i (single-particle states are denoted with 
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small initial letters and many-particle states are denoted 
with capital initial letters.) As we discuss in the follow- 
ing subsection, a crucial observation according to Eq. (||) 
is that there is a degeneracy between states with (n^ — 0; 
m = mo), {ur = 1; m = ttiq — 2), (n^ — 2; m = rriQ — 4), 
etc., with an energy (ttiq -I- l)huj±. States with m having 
the opposite sign as L are excluded, since mixing them 
results in higher-energy |^ Jl^ , |l3|| . 

B. Symmetry- breaking confining potential 

Let us now assume that in addition to the circularly- 
symmetric harmonic potential there is a small symmetry- 
breaking potential AH which has the form 

Ai/ = ^|A/.;i(x?-yn, (3) 

i 

where e <C 1. The eigenvalues of the Hamiltonian H + 
AH are known analytically in cartesian coordinates, and 
they are given by 

^n,,ny = [rixVl + e + ny\/l - e) huj±, (4) 

where and Uy are the number of quanta of oscillation 
along the x and y directions respectively. The corre- 
sponding eigenfunctions are given by the product of the 
eigenstates along the x and y directions with Ux and ny 
quanta, respectively, 

{x/ax)(f>ny{y/ay)- (5) 

Here = (h/MmiY/^ , is the oscillator length along di- 
rection i, with ujx,ujy = uj±^/l ± e. 

Although the non-interacting problem can be solved 
exactly, it is instructive to continue working in cylindri- 
cal polar coordinates, and consider AH as a small per- 
turbation. In these coordinates AH links single-particle 
states which differ by two units of angular momentum, 
m' — m = ±2. Therefore AH links the degenerate states 
with (ur = 0; m — mo), {rir — 1; m — mo — 2), 
{rir = 2; m = mg — 4), etc. Within the subspace of 
these degenerate states the matrix elements of AH are 
all zero, except the ones linking neighboring states, and 
they are of order e. Therefore the problem of diagonal- 
izing H + AH is equivalent to a tight-binding model, 
with a position-dependent hopping integral connecting 
only nearest neighbors. The correction Ae to the single- 
particle energies due to AH is thus of order e, and the 
energy of the single-particle states is 

<^nr,m = <^«r,™ + Ae = [mo + 1 + ©(eiv.m)] huj^, (6) 

where ©(en,-,™) denotes a term of order e which depends 
on Ur and m. The corresponding single-particle states 
are given by 

'/'rv.m = '/'n,.,™ + ^ (^<,m' ) -/-^ > (7) 



where and m' can take all the possible values, with the 
exception ^ rir, and m' ^ m. The results of Eqs. (|) 
and (0) can also be seen by expanding Eqs. (||) and (|) 
respectively in powers of e. 

To be more explicit, the wavefunctions we have to con- 
sider for the lowest-energy state of the system are domi- 
nated by the component having no radial nodes, 

with = 1, 2, . . ., and m' = m — 2, m — 4, . . . 

Equation (|^) implies that the lowest energy of N (non- 
interacting) bosons having L units of angular momentum 
(measured relative to that of the ground state) is 

E{L,N) = [L + NO{eL.N)]riUi_, (9) 

where ©(e^.w) is a term of order e that depends on L 
and N. According to Eq. (||) there is a spreading in the 
energy of the states around the value Lhuj±, which is of 
order NeUuj±. As long as e is small, e ^ 1, the effect of 
this term is not substantial for L being 0{N). 

C. Effect of interactions — Hierarchy of terms 

Turning to the effect of the interactions between the 
atoms, we assume that these are of zero range, 

Vint = ^UoY.6{r^,-r,), (10) 

where Uq = Anh^asc/M is the strength of the effective 
two-body potential, with Osc being the scattering length 
for atom- atom collisions. In the present article we ex- 
amine only the case of effective repulsive interactions be- 
tween the atoms, Osc > 0. 

In the limit that the interactions are weak, the energy 
of the states which are degenerate in the absence of in- 
teractions is spread over a width of order nUo^ where n 
is the atom density and thus the following condition has 
to be satisfied 

nllo < huj±. (11) 

If N is the number of atoms in the trap, and assuming 
that the gas is trapped along the z axis on a length scale 
of order i?^, and in the transverse direction on a length 
scale of order a±, n is of order N/a\Rz; the above con- 
dition can then be written as 

vq < huj^/N, (12) 

where vq ~ Uo/a\Rz. Equation ( [l^ ) implies that the 
scattering length has to scale as 

Osc--^, (13) 
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for the assumption of weak interactions to be valid (for 
fixed u!±.) As a result, the ratio eh,uj±/Nvo is <C 1, and 
therefore 



eh,uj± <C Nvq ^ hLu±. 



(14) 



Under the above conditions the many-body states built 
up by (pQ „ are quasi-degenerate. Therefore one has 
to deal with a problem similar to that of Refs. j^lT^] 
of a symmetric confining potential, with a set of 
(quasi) degenerate states whose degeneracy is lifted by 
the interactions. Under the conditions given by Eq. (14), 



the energy of the lowest state (the so-called yrast state) 
with L units of angular momentum and N particles is 
given by 



£{L,N)^E{L,N) + E,^,iL,N). 



(15) 



Here E{L,N) is given by Eq. (^, and we write it for 
convenience as E{L,N) = (L + aNe)hto±, where a is 
a parameter of order unity that depends on L and N. 
Also £'int(i, N) is the interaction energy due to Vint- For 
£ = 0, £'int(i, N) has been expressed as a power series in 
N and L in Ref. |l|, 

Eint{L, N) = {biN^ + b2NL + bsN + b^L + . . .)«„, (16) 

where bi are parameters of order unity that depend on L 
and N. For e 7^ 0, Eiat{L,N) becomes, to leading order 
in e, 

Eint{L, N) = (l + ce){biN^ + b2NL + 63^^ + biL + .. .)vo, 

(17) 

where c is a a parameter of order unity that depends on 
L and N . The terms proportional to ce result from the 
mixing of states with rtj. ^ (to leading order) , but they 
are not of importance in our problem. 

According to the above equations £(L, N) is given by 
(to leading order) 

£{L, N) = Lnuj±_ + {biN + b2L)Nvo + aNehuj± + 

(18) 

with Lhuj± » N^vo ^ Nehuj±, so the terms of order 
N'^vq in the interaction energy dominate the term in the 
energy due to the symmetry-breaking potential. For a 
given L, the corresponding angular frequency of rotation 
il of the gas is 



1 d£{L,N) 



(19) 



The single-particle density contours that we derive in 
Sec. Ill refer to a frame that rotates with this frequency 
ri, which we assume is also the frequency of rotation Qq 
of the asymmetric potential AH (under equilibrium con- 
ditions.) In Sec. IV we examine the response of the gas 
if the trap stops to rotate, i.e., if flo 0. 



D. Connection with experiment 

Before we proceed, we explore the relevance of our 
study to experiment. As we mentioned in the Introduc- 
tion, in the experiment of Ref. Q , more than 100 vortices 
were created and as argued in Ref. ||l^ , in such a system 
one can easily get to the regime of weak interactions. 

It is even more interesting that in the recent exper- 
iment described in Ref. (24|, it has become possible 
to create Bose-Einstein condensates of ^^Na atoms in 
cigar-shaped traps, where transversely to the long axis 
of the trap the gas is in the lowest harmonic-oscillator 
level. More specifically, the trapping frequencies used 
were u)zI'2,'k — 3.5 Hz and lo_\_/2'k = 360 Hz, which im- 
plies that ~ 11 /im and aj_ « 1 /im. To put as many 
atoms as possible, Gorlitz et al. created a Bose-Einstein 
condensate which was in the Thomas-Fermi regime of 
strong interactions along the long axis, but on the other 
hand it was in the lowest harmonic-oscillator level per- 
pendicular to the long axis. Under these conditions, the 
transverse width of the cloud is on the order of the os- 
cillator length a±, and for the interaction energy to be 
comparable to hu;±, 



(20) 



where is the width of the cloud along the z axis. As 
we mentioned, the cloud is in the Thomas-Fermi limit 
along this direction, and therefore 



From Eqs. ( |20[ ) and ( |2lD we find that 
N, —, 



a±as. 



(21) 



(22) 



which gives N,. ^ 5 x 10*. Indeed, it was confirmed that 
as N varied around this value, the cloud passed from the 
limit of strong interactions to the one of weak interac- 
tions. 

Under the above conditions, the coherence length ^ at 
the center of the cloud is 



Nasca± 



N 



(23) 



For N ^ Nc, a± <C ^, i.e., the coherence length becomes 
larger than the size of the trap and the properties of the 
system under rotation resemble those of superfluid nuclei. 



III. DENSITY CONTOURS 

We investigate now the (normalized) single-particle 
density of the atoms p{r), in the many-body state \ '^l.n) 
with N atoms and L units of angular momentum. 
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p{r) = ^(*L.iv 



(24) 



From the exact solutions that come from diagonahza- 
tion of Mnt in the subspace of degenerate states [^D, 
the single-particle density is always axially symmetric, 
reflecting the symmetry of the confining potential. 

Let us now introduce a mean- field wavcfunction |0,|f2| , 
which is a Fock state expressed as a product of single- 
particle states: 

■^L,N{ri,r2, . . . ,rAr) = V(i"i) X 'ip{r2) . . .'0(rjv)- (25) 

The single-particle states ip{ri) can be expanded on the 
basis of the harmonic-oscillator eigenstates 0o,m(ri) with 
no radial nodes and an angular momentum mh along the 
axis of rotation: 



m— 



Cm</'0,m(l"j)- 



(26) 



In the approach of Refs. [|7|Jl2| the coefRcients c„i in 
Eq. (|2^) are treated as variational parameters with two 
constraints imposed on them: the first is the normaliza- 
tion condition, J2m I'-ml ~ ^' ^^'^ second is that the 
expectation value of the angular momentum be equal to 
/ — L/N, X]m"^l'^mP — L/N. One then calulates the 
interaction energy which is 



1 



(27) 



to leading order in TV, and minimizes it with respect to 
the Cm- As a result of this approach, for repulsive interac- 
tions and within the mean-field approximation a vortex 
array develops as L/N increases ]7|,1^. 

The mean-field wavefunction above has the following 
features or advantages: (i) it gives the exact leading- 
order term of the interaction energy, i.e., the N'^vq [13|, 
(ii) the corresponding single-particle density breaks the 
rotational symmetry, since the single-particle states given 
by Eq. ( p6|) have some axes of symmetry by con- 

struction, and (iii) it is essentially a one-body wavefunc- 
tion, and thus it is easy to work with and to visualize. 
As shown in Ref. there is a freedom in one of the 
phases of the coefficients Cm reflecting the rotational in- 
variance of the confining potential. Therefore even in the 
mean-field approximation one can construct a circularly- 
symmetric state by superimposing the above states cor- 
responding to different angles, and therefore to different 
orientations. 

However, even an infinitesimally small e is enough to 
change p(r) from being symmetric, looking like the one 
given by mean-field, and therefore with a specific ori- 
entation. To see this more clearly, we notice that the 
leading-order N^vq term in the interaction energy given 
by the mean-field approximation is identical to the one 
that comes from the exact solution, as shown in Ref. . 
For any nonzero value of e the term of order Nehw± has 



to become as small as possible, as Eq. ( |I8| ) implies, in 
addition to the leading-order term of the interaction en- 
ergy, i.e., the N'^vq. In order to minimize the Nehuj± 
term, the phase that is free to take any value when e = 
(reflecting the rotational invariance of the confining po- 
tential) has to be chosen appropriately. Therefore for 
e 7^ there will be in general some preferred orientations 
of the cloud corresponding to some (distinct) values of 
the phase which minimize the Ne?iuj± term. These axes 
of symmetry rotate with an angular frequency f2 given 
byEq.(|T|). 

Therefore, according to our analysis, an infinitesi- 
mally asymmetric confining potential makes the circu- 
larly symmetric single-particle density of the cloud un- 
stable against the formation of a vortex array as given 
by the mean-field approximation. This array of vortices 
has been investigated in Refs. [j7|Jl2| (see Fig. 2 in Ref. 
and Figs. 3 and 4 in Ref. pl). 



IV. TIMESCALES 



Related to the above remarks are some timescales 
which we examine now. We consider a problem corre- 
sponding to the experiment of Ref. , where the atoms 
get their angular momentum when they are still in the 
normal phase, above the transition temperature, by ro- 
tating AH with some angular frequency r2(t), which 
slowly approaches the value SIq. The gas is then cooled 
down to the condensed phase, in the presence of the ex- 
ternal rotating drive. The system equilibrates in the ro- 
tating frame, and eventually AH is turned off slowly. Af- 
ter some time the atoms are released from the trap and 
expand, and their density is observed. Since the whole 
process occurs adiabatically, the state of the system is de- 
termined by minimizing the energy in the rotating frame, 
i.e., by minimizing £'{L,N) = £{L,N) — Lhrto. 

Therefore initially, when the gas is in the normal phase, 
it has some angular momentum, which changes during 
the cooling process to the condensed phase. When the 
cooling process is completed the gas is not in equilibrium 
yet, since its angular frequency of rotation Q given by 
Eq. (|l^) is not necessarily equal to the frequency of rota- 
tion of the trap f2o- One therefore has to wait for some 
time for the system to equilibrate, so that fl will become 
equal to r^o- After this time interval the gas is driven by 
the rotating asymmetric trap and it is in a "dynamical" 
steady state. 

The interesting question is how does the gas respond 
if the symmetry-breaking potential stops to rotate. If 
this process occurs adiabatically, i.e., on timescales much 
larger that wj^, the gas will go back to rest. The opposite 
limit, when the trap stops to rotate on a timescale that 
is short as compared to wj^, is non-trivial, however. In 
this case the gas experiences a time-dependent potential. 
This potential induces transitions to states with a dif- 
ferent angular momentum. Because of the symmetry of 
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the perturbation AH we have chosen (i.e., — y^), only 
states which differ by two units of angular momentum 
are linked via AH. The rate F for the system to lose two 
units of angular momentum (due to the torque coming 
from AH that no longer drives the gas) is of order 



r = 



2^|($L^2.jv|Ag|^L.jv)P 

h \A£{L,N)\ 



(28) 



Here \^l-2,n) is the state that couples to the ground 
state I^^l.n) via AH, and A£{L,N) is the energy dif- 
ference between them. The rate T is dominated by the 
transition to the state that is closest in energy to the 
yrast. Since the yrast state consists mostly of single- 
particle states with no radial excitations, the state that is 
closest in energy is the one dominated by single-particle 
states with one unit of radial excitations. The energy 
denominator in Eq. (Bq) is, to leading order. 



\A£{L,N)\^ Nvo, 



(29) 



plus corrections of order Nehuj± . 

To proceed to the calculation of the matrix element 
of Eq. (^), let us denote as the occupancy of the 
state of the harmonic oscillator 0o,m when the system 
is in the yrast state \'^l,n)- Within the mean- field ap- 
proach Njn = iVjcmp, whereas within the exact approach 
one needs to project the yrast state on the basis 0o,m 
and get the corresponding Nm- There is very strong ev- 
idence that the two approaches give the same result for 
Nm |T^.|25fl . The above matrix element is on the order of 



L-2.N 



oo 

-E 

m— 



m+2- 



(30) 



Typically m in the above sum runs over a few states only, 
since states with high values of m are not occupied. Now 
two possibilities need to be considered. Depending on the 
value of the ratio L/N, the sum on the right of Eq. ( |30| ) 
can be either of order iV^/^, or of order N ||7|, p^ , p^ and 
r depends crucially on that. We examine each case sep- 
arately below. 



A. Fast relaxation 

When the yrast state consists of at least two states Nm 
and Nm+2 which are of order N, the matrix element of 
Eq. ^) is of order NehwjL and from Eqs. (||) and (|2|) 
we get that the rate for the system to lose L ~ iV units 
of angular momentum is 



r 

N 



vo 



(31) 



where x{L/N) is a function of L/N , and its value is of 
order unity p6[. 



An important conclusion resulting from Eq. (^Tj) is that 
the timescale for relaxation of angular momentum can get 
long as compared to the period of the trap, provided that 



Wo 



1/2 



(32) 



Since, according to Eq. (12) uq ^ Tilj^_/N for the assump- 
tion of weak interaction to be valid, therefore 



1 



e oc 



7V1/2 



(33) 



for the decay rate to be ^ lo^. Under the conditions 
of Ref. H, £ < 0.01 for = 5 X 10^, which is a 
rather restrictive condition. Realistic symmetric traps 
have anisotropics on the order of 1% [g^i and therefore 
the relaxation rate of angular momentum should be on 
the order of the trap frequency, even for the least asym- 
metric traps. 

The above case of fast relaxation applies to all the yrast 
states with L/N < 2.03, except for the region around the 
single vortex, L/iV « 1 Above the value L/A^ « 2.03 
the system shows a discontiuous transition from a two- 
fold symmetric state to a three- fold symmetric state . 
Also for L/N w 1 the yrast state is dominated by the 
occupancy of a single state, the m — 1 These 
cases are examined in the following subsection. 



B. Slow relaxation 

For L/N > 2.03, the yrast state is three-fold symmet- 
ric; therefore for L/N > 2.03 and for L/N w 1, the 
matrix element of Eq. ( pO| ) is instead given by 



{'Pl-2,n\AH\-^l.n) 
In this case, 



Nehujj 



(34) 



F 

N 



• uj±e 



Nv^' 



(35) 



From Eq. (35) one sees that for the state to be stable. 



£ < 



(Nvo^'^' 



(36) 



which is a condition much easier to satisfy, since although 
£ has to be numerically small, it does not have to scale 
with the number of atoms, A^. 

To summarize the results of this section, the presence 
of AH that rotates with an angular frequency fio results 
in (i) transferring angular momentum, and (ii) favour- 
ing the creation of a vortex array instead of a circularly- 
symmetric ground state. Finally stopping the rotation of 
AH results in a relaxation rate of the angular momentum 
of the gas that depends on the yrast state considered, i.e., 
on the ratio L/N. For L/N < 2.03, except L/N w 1, the 
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rate is low as long as e ~ iV~^/^ and the states are sta- 
ble under such weak perturbations. For L/N k, 1, and 
for L/N > 2.03, the relaxation rate is suppressed by a 
factor of 1/iV, and e has to be numerically small, but 
independent of iV, for these states to be stable. 

The relaxation rate can also be evaluated for higher 
values of L/N . According to Ref. |7j, these states have 
certain symmetry, and as long as this symmetry is not 
two- fold pS], the relaxation rate is given by Eq. (B5|). 



V. CONCLUSIONS 

In conclusion, we have investigated the effect of 
a small asymmetry in the confining potential of a 
weakly-interacting Bose-Einstein condensate, which ro- 
tates bringing the gas into rotation. Three important 
results have come out of this study: (i) the limit of weak 
interactions can be achieved, even under current experi- 
mental conditions, (ii) even in a weakly-asymmetric ro- 
tating trap, in its lowest state the system develops a vor- 
tex array for an effective repulsive interaction between 
the atoms, and (iii) if the gas is no longer driven exter- 
nally, the rate of relaxation of the angular momentum 
depends on the yrast state considered. For the specific 
form of the perturbation AT? (i.e., — y'^), if the yrast 
state does not have some symmetry, or if it has a two-fold 
symmetry, the rate of relaxation of angular momentum is 
large. On the other hand, the unit vortex, or states with 
s-fold symmetry, with s > 3 are expected to be more ro- 
bust against anisotropics of the confining potential, since 
the rate of relaxation of angular momentum is lower by a 
factor of 1/N. These are definite theoretical predictions 
and it is interesting to confirm them experimentally af- 
ter the expected creation of vortices in weakly-interacting 
Bose-Einstein condensates. 
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